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We show that higher-order actions for cosmological perturbations in the multi-field DBI-inflation 
model are obtained by a Lorentz boost from the rest frame of the brane to the frame where the 
brane is moving. We confirm that this simple method provides the same third- and fourth- order 
actions at leading order in slow-roll and in the small sound speed limit as those obtained by the usual 
ADM formalism. As an application, we compute the leading order connected four-point function 
of the primordial curvature perturbation coming from the intrinsic fourth-order contact interaction 
in the multi-field DBI-inflation model. At third order, the interaction Hamiltonian arises purely 
by the boost from the second-order action in the rest frame of the brane. The boost acts on the 
adiabatic and entropy modes in the same way thus there exists a symmetry between the adiabatic 
and entropy modes. But at fourth order this symmetry is broken due to the intrinsic fourth-order 
action in the rest frame and the difference between the Lagrangian and the interaction Hamiltonian. 
Therefore, contrary to the three-point function, the momentum dependence of the purely adiabatic 
component and the components including the entropic contributions are different in the four-point 
function. This suggests that the trispectrum can distinguish the multi-field DBI-inflation model 
from the single field DBI-inflation model. 



I. INTRODUCTION 



Precise measurements of the cosmic microwave background (CMB) anisotropies such as those obtained by the 
WMAP satellite [l[ provide valuable information on the very early universe. Any theoretical model that attempts to 
explain the evolution of the universe before the big bang nucleosynthesis will also have to explain the observed CMB 
anisotropies. Even though these anisotropies are almost Gaussian, a small amount of non-Gaussianity is still allowed 
by the data USUI. The information contained in this non-Gaussian component will contribute to a huge advance in 
our understanding of the very early universe. For example, the simplest slow-roll single field inflation models predict 
that the non-Gaussianity of the fluctuations should be very difficult to be detected Q, even in future experiments 
such as PLANCK [6( . If we detect large non-Gaussianity, this means that the simplest model of slow- roll single field 
inflation would be rejected. 

Recently, theoretical models which can produce sizeable non-Gaussianity has been extensively studied by many 
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relax one (or more) of the following standard conditions: single field, slow-roll, canonical kinetic term and standard 
Bunch-Davies vacuum. In the present work, we will consider the case of multiple fields with non-canonical kinetic 
terms. For example, in i^-inflation models where the kinetic term of the infiaton field is generic, the sound speed of 
the perturbations can be much smaller than 1 [62 . 63], which leads to large non-Gaussianity. 

Among the models, the Dirac-Born-Infeld (DBI) inflation, motivated by string theory, can also realize large non- 
Gaussianity [HI, [I?], 0, HH, [ill [6?| • In this model, the infiaton is identified with the position of a moving D3 brane 
whose dynamics is described by the DBI action. Again, due to the non-trivial form of the kinetic term, the sound 
speed can be smaller than 1 and the non-Gaussianity becomes large. However, recently it has been pointed out that 
DBI-inflation driven by a mobile D3 brane with large non-Gaussianity might contradict the current WMAP data. For 
current and stringent observational constraints and consequences of DBI-inflation see (gH, H^, [7(1 [Zll, [z3, [zD, [z3) [zH 

One way to avoid this constraint is to consider multi-field DBI models (36| . Since the position of the brane in each 
compact direction is described by a scalar field, DBI-inflation is naturally a multi-field inflationary model [78j]. As 
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first pointed out by [79j], in multi-field inflation models, the curvature perturbation is modified on large scales due to 
the entropy perturbation. Even though there are some works considering multi-field inflationary models with kinetic 
terms depending on X — —Gijd )Jl 4> I d fJ '(j) J /2, where ip 1 are the scalar fields (/ = 1, 2, ...) and Gu is the metric in the 
field space as in the case of K-inflation (4f| [8(| HH , the consistent analysis for the entropy modes in the multi-field 
DBI- inflation model has started only very recently [35|, [H, [U HH, [57j • 

In [H, H(| [37], the three-point function in the small sound speed limit and at leading order in the slow- 
roll expansion was obtained and it was shown that in addition to the purely adiabatic three-point function, there 
exists a mixed component (Q (T (ki)(3, j (k2)(3 s (k3)) where Q a and Q s are the adiabatic and the entropy perturbations, 
respectively. Since the momentum dependence of the three-point function from the adiabatic modes was shown to be 
identical with the mixed component, the shape of the bispectrum of the curvature perturbations remains the same as 
in the single-field case, while the amplitude is affected by the entropy perturbation. 

The previous works on the non-Gaussianity in the multi-field DBI-inflation model are limited to the bispectrum, 
with the exception of [55[ , where the authors compute the leading order trispectrum based on the assumption that it 
is mainly from some limited terms of the entropy perturbations. It is expected that the cosmic microwave background 
(CMB) trispectrum also provides strong constraints on early universe models. At the moment, the constraints are 
rather weak given by |tjvl| < 10 s [82l |83||. where tml denotes the size of the trispectrum. However, PLANCK will 
tighten the constraints significantly reaching \tnl\ ~ 560 84]. It is also worth noting that the analysis in the single 
field DBI-inflation model shows that the trispectrum is enhanced in the small sound speed limit as tml ~ 1/c f 
[111, [HI, [86[ . As in the bispectrum case, the constraints depend on the shape of the wave vectors' configuration pl\ - 
Therefore, it is important to calculate the shape dependence of the trispectrum in the multi-field DBI-inflation model. 
For the details of the observations of the CMB trispectrum, see (88l. 183. |90]| . 

In this paper, we calculate the four-point function of the primordial curvature perturbation coming from the intrinsic 
fourth-order contact interaction in the multi-field DBI-inflation model and see whether the momentum dependence 
of the four-point function is useful to discriminate the multi- field DBI-inflation model or not. In order to obtain the 
third- and fourth- order actions for cosmological perturbations, we propose a simple and intuitive method based on 
a Lorentz boost, making use of the special property of the DBI action. This does not only provide the fourth-order 
action easily, but also explains why the momentum dependence of the three-point function from adiabatic modes is 
identical with the one from the mixed component. 

The structure of this paper is as follows. In section [III we describe our model and define perturbations in the flat 
gauge. We decompose the perturbations into the adiabatic and entropy directions. In section IIIII after obtaining 
the fourth-order action at leading order in slow-roll and in the small sound speed limit in terms of the decomposed 
fields based on the ADM formalism [f| [l?], Oil [9l| , we calculate the fourth-order interaction Hamiltonian. Then, the 
four-point functions coming from the intrinsic fourth-order contact interaction are derived in section IIVI In section 
fVl in order to develop understanding of the results, we derive the same fourth-order action based on a Lorentz boost 
which relates the brane-rest frame and the brane- moving frame. Section [VTl is devoted to the conclusion. 



II. THE MODEL 



We start with the multi-field DBI-inflation model described by the following action [92 



S 



P(x,4> 1 ) = - 



R + 2P(X,(t) 1 ) 



1 - 2f(tp)X - 1 - Vitf) 



(1) 



where we have set 8nG = 1, R is the Ricci scalar, <// are the scalar fields (I — 1,2, ...,N), /(0 7 ) and V{(f) 1 ) are 
functions of the scalar fields determined by string theory configurations and X is defined in terms of the determinant 
V = det(<5£ + fGud^tfdvtf) as X = (1 - V)/(2f). Here Gu is the metric in the field space. We assume that P is 
a well behaved function of (j) 1 and X. It is also shown that X is related to the kinetic terms of the scalar fields as 

EE El 



X = G U X IJ - 2fx} I Xf ] 



ifxfx/Xp - 8fX r [I Xj J X^X 



L] 

L > 



X 



IJ 



x/ 



r< yKJ 



(2) 
(3) 



where the brackets denote antisymmetrization. It is worth noting that even though X = X ( = GijX ij ) in the 
homogeneous background, this does not hold if we take into account the inhomogeneous components. 
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In the background, we are interested in flat, homogeneous and isotropic Friedman-Robertson- Walker universes 
described by the line element 

ds 2 = -dt 2 + a 2 (t)5 i] dx t dx j , (4) 
where a(t) is the scale factor. The Friedman equation and the continuity equation read 

3H 2 = E , (5) 



E a = -3H \E + PqJ , (6) 

where the Hubble rate is H — a/ a, dot denotes derivative with respect to cosmic time t, Eq is the total energy of the 
fields which is given by 

E = 2X{ ) J Po,x" -Po, (7) 

and the subscript denotes that the quantity is evaluated in the background. 

For the later convenience, we introduce the following parameter that characterizes the motion of the brane in this 
background and serves as a Lorentz factor 

1 



7 (<^,X ) = -==, with v = ^2f X Q , (8) 

where we have used the fact that Xq — Xo. 

For this model the speed of propagation of the scalar perturbations ("speed of sound"), c s , is given by 

W fi ) (9) 



Pf + 2XP xx 



n 



It can be shown that c s is the inverse of 7. 

Since we are interested in the inflationary background, we assume the form of /(</> J ) and V{<p I ) are chosen so that 
the inflation is realized at least for 60 e-foldings. In order to characterize this background, we define the slow-variation 
parameters, analogues of the slow-roll parameters, as 

H X Q e c s 

e = -Jp = H 2 V s ^ v = 7h> x = ^h- (10) 

We should note that these slow-variation parameters are more general than the usual slow-roll parameters and that 
the smallness of these parameters does not imply that the field in rolling slowly. We assume that the rate of change 
of the speed of sound is small (as described by \) but c s is otherwise free to change between zero and one. 

We shall consider perturbations on this background. We decompose the scalar field ip 1 into the background value 
4>q and the perturbation Q 1 in the flat gauge as, 

<f> I (x,t)=4> I Q {t)+Q I {x,t). (11) 

Furthermore, as was done in [93|, we decompose the perturbations into instantaneous adiabatic and entropy per- 
turbations, where the adiabatic direction corresponds to the direction of the background fields' evolution while the 
entropy directions are orthogonal to this. We introduce an orthogonal basis e J n , with n = 1, 2, N, in the field space 
so that the orthonormality conditions are given by [35[ 

/ 1 1 - c l 
e„e m / = — S mn -SmiSni , (12) 

Cs C s 

where the adiabatic basis is defined as 



— 4>i. (13) 
2X V ' 

Notice that the length of the basis vector e{ is c s and that of the other basis vectors is l/c s - If we consider the 
two-field case (I = 1,2), the field perturbations are decomposed into the adiabatic field Q a and the entropy field Q s 
as 

Q 1 = Q*e[ + Q s ei . (14) 

Hereafter, for simplicity, we will concentrate on the two-field case although the extension to more fields is straightfor- 
ward. 
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III. FOURTH-ORDER ACTION AND HAMILTONIAN AT LEADING ORDER 



We first calculate the fourth-order action for Q n , where the subscript n denotes either a or s. Since we are interested 
in the leading order shape of the trispectrum, we keep only the leading order in the slow-roll approximation, where the 
values of the slow-variation parameters defined by Eq. (fTU|) are assumed to be small. We also assume c s <C 1 because 
otherwise the trispectrum is not observable in the future experiments. Using these approximations and following the 
ADM formalism 0, [13, [H, Hl| , the action up to fourth order can be calculated as 



S, 



( m t 
(2) 



s. 



(main) 
(3) 



, or 

dtcPx^r 

d 



= - / dtd A x- 



S, 



(main) 



(4) 



1 

16 



Ql + Ql - ^ ( d t Q a d l Q a + d t Q s d l Q t 



Ql + Q a Q 2 s + % ( (diQ s d l Q s - diQ^Qa) Q a - 2 (d.Q^Qs) Q t 



„3 



= — / dx A dt 



clX 



2c 



(15) 
(16) 



+Q 2 AQad i Q„ + Q 2 AQ s d l Q s 
m {d l Q a d l Q a ) 2 -2{d i Q a d i Q„) (djQsd'Qs) +4(d l Q a d l Q s ) 2 + (d t Q s d l Q s y 



(17) 



In general, the interaction Hamiltonian is not just the opposite sign of the interaction part of Lagrangian as is 
explained by the following calculation. Following |28| , we define the Lagrangian density as 



c = /< o) a?+/r^+/ 2) +5i o) ^^^ 



(18) 



where a m with m = 1,2 denotes Q a and Q s respectively. The /'s, j's, g's and /i's are all functions of a m (t,x), its 
spatial derivative dia m and time t. They come from 2nd, 3rd and 4th order action, respectively. For these functions, 
the superscripts (0), (1), (2), (3), (4), denote the order in the perturbations a m of these functions. The subscripts a, 
b, c, just label the functions. 

The momentum densities are given by 



Making use of these relations, a m are expressed in terms of 7r m up to third order as 



Ql 



~1 



2/i 



(0) 



c^l+c^n 2 



a 2 



7T 2 



2/, 



(0) 



di 



4V 



7T2 



4 3) -l. 



(19) 
(20) 

(21) 



where the different c's and d's functions are given by 



,(2) 



„(3) 



dl 3) = - 





(2 
9 a 


(2) 


8(/i 0) )3 2/i°> 


2' °b 




hg a %^ 


■ 9b a a 


2/i 0) 


V /i 0) 


+ f (0) 
Jb 




( 9 f ] d a 2) 


o (0) c (2) 

9b °a 


2/ b (0) 




r AO) 

Jb 




W 




2/ b (0) 


\ fl 0) 


f 2(/D 3 



(0) 

9 b 



fi 0) (/ 6 (0) ) 2 



,(3) 



2/i' 



(0) 



,(()) 



,(2) 



4/i 0) (/ b (0) ) 2 , 

(0)x 2 f(0) I 



. ,(0) »(Q) 
4/o 7fc ^1^2 



,(2) 



,1 j-(O) j-(0) 

4/a 7 6 ^1^2 



,(2) 



2(/ fc (0) ) 2 - 2 2 



f (o) 



(0) 



2(/A' 



(0)> 



■9b 



(0) 



(2) 



2(/i 0) ) 2 7T 2 



36 



(2) 



,(22) 



,(23) 



(24) 
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where again for these functions the superscripts (0), (1), (2), (3), denote the order in the perturbations a m of these 
functions and the subscripts a, b, just label the functions. 

The Hamiltonian density is obtained by plugging the previous expressions into 

H = 7Tidi + 7r 2 d 2 — £ • (25) 

We then separate Ti into a kinematic Hamiltonian density Ti$ which is given by 

n °-TiF + W'- l{2K <26) 

and an interaction Hamiltonian density Tt mt . To use the interaction picture formalism [94| . a m and 7r m in this 
interaction Hamiltonian should be replaced by their interaction picture counterparts and ir^, which satisfy the 
free equation of motion determined by 7i - They also satisfy the usual commutation relations 

[c£ (t, x) , ni (t, y )] = iS mn S 3 (x - y) . (27) 

Expressing in H mt in terms of d^ using 

j dH 7rf . <9W , 9sn 

we finally get the third-order and fourth-order interaction Hamiltonian densities as (omitting the label "I" in the 
variables in Ti mt from now on) 

= ~ - - flP'ia - J (3) , (29) 



UT{ = Hg^l (0) 3flW + OT _ h m ^ + f MT _ ,(0) d4 

(4) " l 4/W ° J 1 I 2/W + /<°> " b J ^ + I 4/^ " C J " 2 



3ga ga _ U? + gb gb - h 2) I dido + I gb ga - hP> d 2 + (ga } + {9b ' - 7 < 4 ) 



(30) 



As in the single-field case [28| , while the cubic part of 7i mt is the opposite sign of the cubic C mt , this is generally not 
true at fourth order. As we will see, the extra terms contribute to the leading order results for c s « 1. By applying 
Eq. (f3"0")) to the two-field DBI- inflation model given by Eqs. (JT5J) , (|16p and (|T7|) . we can obtain the fourth-order 
interaction Hamiltonian density as 

m = + ^& 2 s + si? & Q ° diQ ^ « + db? « ■ (31) 

IV. THE LEADING ORDER IN SLOW-ROLL FOUR-POINT FUNCTION 

In this section, we derive the connected four-point functions of the adiabatic and entropy fields coming from the 
intrinsic fourth-order contact interactions at leading order in the slow-roll expansion and in the small sound speed 
limit. The perturbations are promoted to quantum operators as 

Qn(r, x) = ^ I rf 3 kQ„(r, k)e lk x , (32) 

where r denotes conformal time and 

Q„(t, k) = u„(r, k)a„(k) + <(t, -kjaj, (-k) , (33) 

and a„(k) and a^(— k) are the annihilation and creation operators, respectively, that satisfy the usual commutation 
relations: 

[a n (k 1 ),a+ 1 (k 2 )] = (2Tr) 3 S^(k 1 -k 2 )J nm , [a„(k 1 ),a m (k 2 )] = [ a t (ki), at,(k 2 )] =0. (34) 
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At leading order, the solutions for the mode functions are given by 



u n (r,k) =A n ^(l + ikc s T)e- ikc > T . (35) 



The two-point correlation functions are then obtained as 



2n 2 

(0|Q n (r = 0,k 1 )Q m (r = 0,k 2 )|0) = (2 7 r) 3 ^(k 1 + k 2 )V Qn - TT S nm , (36) 

K 1 



where the power spectra VQ n are defined as 



and they should be evaluated at the time of the sound horizon crossing c s ^k\ = a»i/*. 

In terms of these quantum operators, the connected four-point correlation function coming from the contact inter- 
action in the interaction picture formalism is given by (f| l94j 

t 



(n|Q m (t,k 1 )Q n (t,k a )Q J ,(t,k 8 )Q,(t,k 4 )|n) = -i I dt(o\ [Q m {tM)Qn(tM)Q P (tM)Q q (tM),H mt (t)] |o), 

(38) 



to 



where to is some early time during inflation when the fields' vacuum fluctuations are deep inside the horizon and t is 
some time after the horizon exit. |f2) is the interacting vacuum which is different from the free theory vacuum |0). If 
one uses conformal time, it is a good approximation to perform the integration from — oo to because r w — (aH) . 
The interaction Hamiltonian H mt is given by H mt = J d 3 xTi mt . 

The purely adiabatic, purely entropic and mixed components are given by 



H* 1 



(fi|g ff (0,k 1 )g ff (0,k 2 )Q CT (0,k 3 )Q CT (0,k 4 )|fi) = ^W^k;) 3 (-36^) , (39) 

TT 4 k 2 

A = %^, * = (40) 



K 5 

i=l 



4 ^8 



(fi|g s (0,k 1 )g s (0,k 2 )Q s (0,k 3 )Q s (0,k 4 )|fi) = (^fS^^kO ^ n4 ^ 3 (~^aj , (41) 



fcffcf(k 3 -k 4 ) / 3(fc 3 + fc 4 ) 12fc 3 fc 4 \ 
A2 = «5 ( v 1 + ^^ + ^^J +perm - 



(42) 



' -f / f 



(n|Q ff (o, ki)Q g (o, k 2 )g s (o, k 3 )g s (o, k 4 )[o) = (2^) 3 ^ 3 >(]T k t ) 2X ^ n4 (Mi + j , (43) 

_ fc^ 2 2 (k 3 -k 4 ) / 3(fc 3 + fc 4 ) 12fc 3 fc 4 \ , , 

3 JT 3 V K K 2 J ' 1 j 

where "perm." in Eq. (|42[) denotes the 23 permutations of the four-momenta. The purely adiabatic component agrees 
with the result of the single- field DBI-inflation model [281 ] . 

We need to relate the four-point functions of the scalar fields to the four-point function of the comoving curvature 
perturbation 1Z which is closely related to the observable quantity. As in [35|], 1Z and Q a are related as 

n= ^r Q - (45) 

and it is convenient to define the entropy perturbation S as 

so that the power spectra are Vs„ — Viz t , where the subscript * means that the quantity should be evaluated at the 
sound horizon crossing. 
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In this work, we ignore the possibility that the entropy perturbation during inflation can lead to a primordial entropy 
perturbation that could be observable in the CMB. But we shall consider the effect of the entropy perturbation on the 
final curvature perturbation. Following the analysis of [sB|, we describe the conversion of the entropy perturbation 
into the curvature perturbation by a transfer coefficient T-jis- Then the final curvature perturbation is expressed in 
terms of the adiabatic and entropy field perturbations as 

TZ = A a Qo* + A S Q S * , A a = (^=) > A s = T ns (^=J • (47) 

Hence the connected four-point function of 1Z at leading order is given by 

(^(k 1 )^(k 2 )^(k 3 )^(k 4 )> - ^(<3 CT (k 1 )Q CT (k 2 )Q <T (k 3 )Q CT (k 4 )) +^^((Q CT (k 1 )Q <T (k 2 )Q s (k 3 )Q s (k 4 ))+perm.) 

+^<Q s (k 1 )Q s (k 2 )<2 s (k 3 )Q s (k 4 )) , (48) 

where "perm" denotes five permutations of the four-momenta. This constitutes one of the main results of this work. 

Since the mixed component and the purely entropic component have different momentum dependence from the 
purely adiabatic component, the momentum dependence of the resultant four-point function of the curvature pertur- 
bation is different from that in the single-field DBI-inflation model [28|. The only effect for the bispectrum due to 
the presence of the multiple fields is a change in its amplitude with respect to the single field case, however for the 
trispectrum the presence of the multiple fields affects also the shape dependence. So in principle, the trispectrum can 
be used to distinguish between the multi-field and the single- field DBI-inflation models. 



V. LORENTZ BOOST 

In this section, we will present an alternative and simpler method to obtain the leading order 1 action for the 
perturbations. 

It is well known that the DBI action (TT]) describes the motion of a brane in a higher dimensional spacetime. For 
simplicity, let us take g^ v as the Minkowski metric 2 and / is constant. Then in the frame where the brane is at rest in 
the background (brane- rest frame), X will be a small quantity because it is written in terms of perturbations of the 
brane positions. Then we can expand the Lagrangian fl}, by ignoring the potential terms which do not contribute at 
the leading order as 

P = X - f -X 2 + fXjX'j + O(Xfj), (49) 

where X and Xu are written in the coordinates of this new frame that we will denote by (t, x, a, s). a is the 
coordinate along the direction of the motion of the brane in the background (adiabatic direction) and s parameterizes 
the orthogonal direction (entropy direction). In the brane- rest frame, the DBI action for the two-field model can be 
written in a fairly simple form as 



<£dx,(50) 



where we have ignored higher-order terms. 

We are interested in the behaviour of the perturbations in the frame corresponding to the set-up shown in the 
previous sections and for this purpose, it is necessary to know the Lagrangian of the perturbations in such a frame. 
From Eq. ((5]), this frame is the one where the brane is moving with the velocity v = ^/2fXo in the background 
(brane-moving frame). If we label the coordinates of the brane- moving frame as (t,x, er, s), the coordinate variables 
in these two frames are related by a Lorentz transformation 

t = l(t + v^~fa), V / 7^ = 7(v / 7^ + ^~), (51) 



1 In fact, with this method one can also obtain the sub-leading terms in c 2 s . 

2 In our case g^ v will be an inflating FRW metric but this is conformally Minkowski and the final result will be the same as for the 
Minkowski case, up to powers of the scale factor. 
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and its inverse transformation 



t = 7 



-vt) , 



(52) 



with x and s unchanged. From Eq. (|5ip . in terms of the background value <jq, v can be expressed as v 2 = f&Q, which 
means <jq — V%Xa- We are interested in the small sound speed limit or equivalently when the brane is relativistic, 
i.e. v ~ 1. 

In order to see the behaviour of the perturbations in the brane-moving frame, it is convenient to introduce a new 
variable as 



6a(t, x) = ait, x) 



vt a(t,~x.) 

7f = ~T : 



(53) 



which is nothing but the adiabatic perturbation in the brane-moving frame, since vt/y/f can be interpreted as ctq in 
the case &q is constant. For the entropy perturbation in the brane-moving frame, we can continue to use s, since s 
is invariant under the Lorentz transformation. It is worth noting that these 5a and s are related with Q a and Q s 
introduced in Eq. (|14[) as 



5a 



t'sQa-i 8 — , — s- 



(54) 



since Q a and Q s are defined in terms of the basis satisfying (fl~2)) and (fl~3|) . 

For the Lagrangian in the brane-moving frame it is necessary to know the transformation law of not only the 
coordinate variables, but also its derivatives. Since the time derivative of a quantity q in the brane-moving frame is 



Of 



M)-|«^-i«-)(i-^)-r-'|«-), 



where we define 7 as 



7 = fl.-vy/f-yS&j ^ 1 (l + v^] 1 2 5a + v 2 f 1 i (5a) 2 ) +o(^5a) 3 ^j : 



(55) 



(56) 



we can express the time derivative of a quantity q in the brane-rest frame by the one in the brane-moving frame as 

(57) 



dq a \ ~ d Q t . . 
9f(t,x)= 7 ^(t,x). 



Similarly, the spatial gradient of a quantity q in the brane-rest frame can be related with the one in the brane-moving 
frame as 



Vq(i, x) = Vq(t, x) + vy/f^V(6a) 



(58) 



One can now use Eqs. (|5D|) . (152p . (fB"4"l) . (f5T|) and |58p to obtain the actions up to fourth-order in the brane-moving 
frame as 



~ / — 

" ! 2c 2 



Qa + - 4 {diQad l Q a + d l Q s d i Q s 



dtdn, 



(59) 



5(3) 



2^2X^1 



2 



dtdx., 



(60) 



C(4) 



16 C 5X 

2c 2 



5Q4 + QQlQl + Ql 



4 

- J {faQ^Qv) 2 - 2 (diQ^Qv) (d.Qs&Qs) + 4 (d,;Q CT d l Q s ) 2 + (d t Q s d*Q s ) 2 } dtdx, (61) 
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where we introduced the scale factor dependence, we used dt = dt/j, v ~ 1 and / ~ 1/oq ~ 1/(2Xq). Although cubic 
interactions are absent in the brane-rest frame (see Eq. (|50[) ). they are induced by the boost in the brane-moving 
frame (see Eq. (|60| ). At fourth order, quartic interactions are already present in the brane-rest frame but additional 
terms are induced in the brane-moving frame. It is worth noting that Eqs. (|59[) . ([SO]) and (fBTj) agree with Eqs. (|15[) . 
([16]) and (fl7|) derived by the usual ADM formalism. 

Now we are in the position to discuss the momentum dependence of the bispectrum and trispectrum. The second- 
order action in the rest frame of the brane is symmetric under the exchange between the adiabatic and entropy modes. 
At quadratic order, because of the change of the time coordinate, dt « jdt, the sound speed in the brane-moving 
frame deviates from the speed of light. Since the shift of the sound speed is solely due to this coordinate change, 
both the adiabatic and entropy modes should have a common sound speed [37j . The third-order action in the brane- 
moving frame, which is generated by the boost transformation, also originates from the second-order action in the 
brane rest frame. As is seen from Eqs. (57) and (58), the boost acts on both modes in the same way. Therefore the 
boost transformation preserves the symmetry between the adiabatic and entropy modes in the sense that the mixed 
component which contains Q s reduces to the purely adiabatic component if all Q s are replaced with Q a . At third 
order, the interaction Hamiltonian is just the minus of the third-order interaction Lagrangian. Thus, at this order, 
the three-point function from the mixed component has the same momentum dependence as in the purely adiabatic 
case. Thus, adding the contributions from the entropy modes does not change the momentum dependence. 

At fourth order, there are two contributions to the action in the brane-moving frame. One is the fourth-order action 
that arises by the boost from the second-order action in the brane-rest frame: 



S, 



(4) 
boost 



4c5A 



dtdx. 



(62) 



As in the case of the third-order action, the boost preserves the symmetry between the adiabatic and entropy modes. 
Thus, the mixed component becomes the same as the pure adiabatic component by replacing Q s with Q a and the mixed 
component would give the same shape dependence of the trispectrum as the pure adiabatic component. However, at 
fourth order, there exists the intrinsic fourth-order action in the brane-rest frame which contains pure adiabatic, pure 
entropy and mixed components. If one adds these contributions to Eq. (|62|) . the symmetry between the adiabatic 
mode and entropy mode is broken. Moreover, at fourth order, the interaction Hamiltonian is not simply the opposite 
sign of the fourth-order Lagrangian and the additional contributions in the interaction Hamiltonian also do not have 
the symmetry between the adiabatic mode and entropy mode. Thus the momentum dependence of the trispectrum 
becomes different for adiabatic, entropy and mixed components. 

Therefore the third order is a special case where the shape of the bispectrum is not changed by the entropy modes. 
This is because, at third order, the interaction Hamiltonian only arises by the boost which preserves the symmetry 
between the adiabatic and the entropy modes. At higher orders, we expect that the momentum dependence of the 
n-point function induced by the entropy modes is different from the pure adiabatic contribution. This would be crucial 
to distinguish between single field and multi-field models by the shape dependence of the n-point functions. 



VI. CONCLUSIONS 



In the multi-field DBI-mfLation model, it had been shown that the sound speeds for the adiabatic and entropy 
perturbations are the same [36j . It was also shown that the momentum dependence of the three-point function of the 
final curvature perturbation remains the same as in the single-field case, that is, the components of the three-point 
function including the entropy perturbations only change the amplitude of the three-point function from the purely 
adiabatic component. This is because there exists a symmetry under the exchange between the adiabatic and entropy 
modes in the second- and third-order actions for Q n at leading order in slow-roll and in the small sound speed limit. 

In this paper, as a natural extension of these works, we studied the non-Gaussianity from the contact interaction 
trispectrum in the multi-field DBI-inflation model. We first derived the fourth-order action for the perturbations 
based on the usual ADM formalism fEq. (|17|) V It is written in terms of the adiabatic and entropy perturbations in 
the small sound speed limit and at leading order in the slow-roll expansion. 

After deriving the relation between the interaction Lagrangian and the interaction Hamiltonian, which can be 
applied to a fairly general two-field model, we obtained the fourth-order interaction Hamiltonian (Eq. (|31|1 ). It is 
worth noting that while the cubic part of the interaction Hamiltonian is the opposite sign of the Lagrangian density, 
this is generally not true at higher orders [28| . 

Using this interaction Hamiltonian, we derived the connected four-point function coming from the intrinsic fourth- 
order contact interaction, in the small sound speed limit and at leading order in the slow-roll expansion. In these 
approximations, in addition to the purely adiabatic four-point function (Q CT (ki)(3 (T (k2)(3 -(k3)(2 cr (k4)), there exists a 
purely entropic component (Qs(ki)Qs(k2)Qs(k3)(5 s (k4)) and a mixed component (<5o-(ki)Qcr(k2)(5 s (k3)(5 s (k4)). It 
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was shown that the purely entropic and the mixed components have different momentum dependence from the purely 
adiabatic component (Eqs. (|4T) and (|4U)) ). Because of this it was shown that the momentum dependence of 

the four-point function of the comoving curvature perturbations is affected by the entropy modes (Eq. ([48]) ). and the 
shape is different from the single-field case. In contrast to the shape of the bispectrum, which does not distinguish the 
multi-field DBI-inflation model from the single-field DBI-inflation model, the CMB trispectrum can provide a useful 
discriminator for the multi-field DBI-inflation model. 

We also derived the fourth-order action for the perturbations by an alternative and simpler method. Since the DBI 
action describes the motion of a brane in a higher dimensional spacetime, the action for the perturbations can be 
obtained by a Lorentz boost from the frame where the brane is at rest in the background to the frame where the brane 
moves at the velocity v with v 2 = 2f X = /oCTq. In the small sound speed limit (v — > 1), the actions up to fourth 
order in the frame where the brane is moving were calculated (Eqs. (JHU), ([SO)) and (ffJTj) ) and it was shown that they 
coincide with the ones obtained by the usual method (Eqs. (fT5|) , (fTI))) and (fT7| ). From this derivation, we found that, 
at third order, the interaction Hamiltonian arises purely by the boost and it has the symmetry under the exchange 
between the adiabatic mode and entropy mode. This is because the boost does not distinguish between them. At 
fourth order, there exists the intrinsic fourth-order action in the rest frame of the brane. This breaks the symmetry. 
In addition, the interaction Hamiltonian is not just the opposite sign of the Lagrangian. The additional terms in 
the interaction Hamiltonian also break the symmetry. Hence the trispectrum coming from the entropy modes has a 
different shape dependence from the one coming from the adiabatic modes. 

In order to calculate the effect of the entropy perturbation on the curvature perturbation, we need to specify a 
model that describes how the entropy perturbation is converted to the curvature perturbation. In this paper, we 
modeled this transfer by a transfer function Tus ■ It would be interesting to study this mixing in specific string theory 
motivated models. 

In this paper, we considered the trispectrum coming from the intrinsic fourth-order contact interaction. However, 
as it was shown recently in HH, [86j , there are other important contributions for the t risp ectrum coming from the 
interactions at a distance such as the exchange of scalar particles. In particular, Ref. [5a] calculated a part of the 
trispectrum coming from the entropy modes, i.e. (TZ(k 1 )TZ(k2)TZ(ks)TZ(k4,)) cx T^ iS (Q s (ki)Q s (k2)Q s (k 3 )(3 s (k4)), 
when an adiabatic scalar particle is exchanged. We will present the full leading order trispectrum in the multi-field 
DBI-inflation model in a separate publication [9t| . 
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